Hydrodynamical description of a hadron-quark first-order phase transition. 
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Solutions of hydrodynamical equations are presented for the equation of state of the Var der 
Waals type allowing for the first order phase transition. Attention is focused on description of the 
hadron-quark phase transition in heavy ion collisions. It is shown that fluctuations dissolve and 
grow as if the fluid is effectively very viscous. Even in spinodal region germs are growing slowly due 
to viscosity and critical slowing down. This prevents enhancement of fluctuations in the near-critical 
region, which is frequently considered as a signal of the critical point in heavy ion collisions. 
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There are many phenomena, where first-order phase 
transitions occur between phases with different densities. 
Description of such phenomena should be similar to that 
for the gas-liquid phase transition. Thereby it is worth- 
while to find corresponding solutions of hydrodynami- 
cal equations. Though some simplified analytical [l], 3 
and fragmentary two-dimensional numerical solutions 3J 
have been found, many problems remain unsolved. In nu- 
clear physics different first-order phase transitions (e.g., 
to pion, kaon condensates and to the quark state) may 
occur in neutron stars 0, 0| and in heavy ion collisions 
0, [1] . At low energies gas-liquid transition occurs @ . It 
is also expected that at finite baryon density the hadron - 
quark gluon plasma (QGP) phase transition, which might 
manifest itself in violent nucleus- nucleus collisions, is of 
the first-order [6(. The hydrodynamical approach is ef- 
ficient for description of heavy-ion collisions in a broad 
energy range (e.g. see [SHIl])- 

In this letter the dynamics of a first-order phase transi- 
tion is described by equations of non-ideal non-relativistic 
hydrodynamics: the Navier-Stokes equation, the continu- 
ity equation, and general equation for the heat transport. 
We solve these equations numerically in two spatial di- 
mensions, d = 2, and analytically for arbitrary d in the 
vicinity of the critical point. Then we perform estima- 
tions for the case of the hadron - QGP transition. 

The best known example to illustrate principal features 
of a first order phase transition is the Van der Waals fluid. 
The pressure is given by Pyw [n, T] — nT / (1 — bn) — n 2 a, 
where T is the temperature, n is the density of a con- 
serving charge (e.g., the baryon charge), parameter a 
governs the strength of a mean field attraction and b 
controls a short-range repulsion. In practice we use 
a modified Van der Waals (mVW) equation of state 
(EoS): P[n,T] = f(T)P vw [n,T], where the function 
f(T) is chosen so that §g? | (n er ,T or ) = at the critical 
density n cr and the critical temperature T cr . We use 
f(T) ~ C(ST) [1 - 2(6T)} 2 , ST = (T - T cr )/T cr , and 
the pre-factor C(ST) = [1 + 4((5T) 2 ]~ 1 is chosen to re- 
produce ideal gas EoS for sufficiently low n and high 
T, This modification allows us to parameterize the EoS 



with two minima in the free energy, being convenient 
for analytical treatment of the problem near the criti- 
cal point (n cr ,T cr ). We expand the quantities entering 
EoS and equations of hydrodynamics near a reference 
point (p r ,T T ) chosen somewhere in the vicinity of the 
critical point on the plane P{p, T), where p — mn is the 
mass density, m is the mass of the constituent. Assuming 
smallness of the velocity u(f, r) of the germ we linearize 
hydrodynamical equations in u, density Sp = p — p r and 
temperature ST = T — T r . Applying then operator "div" 
to the Navier-Stokes equation and taking z = divu from 
the continuity equation we obtain [l], [2J : 
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d = 2(d- I)/ d. The pressure SP = P - P[p r ,T r ] = 
P ^ F S(S P )'^ \ T * s GX P rcssc d through the Helmholtz free 
energy F for slightly inhomogeneous configurations; r\ x 
and Cr are the first (shear) and second (bulk) viscosities; 
.. + d 2 . 

Note that thus derived Eq. differs from the phe- 
nomenological Landau equation for the nonconserving or- 
der parameter d t 4> = — l{&F / °"4>)i 7 = const, and from 
equations used for the description of the dynamics of 
first-order phase transitions in heavy ion collisions Q 
and in relativistic astrophysical problems The dif- 
ference with the Landau equation disappears, if one sets 
zero the square bracketed term in the r.h.s. of Eq. |T]). 
From the first glance, such a procedure is legitimate, if 
space-time gradients are small. However for a germ, be- 
ing prepared in a fluctuation at t — with a distribution 
Spit = 0, r) = 5p(0.r), the condition dSp(t, r)/dt\ t= Q ~ 
should also be fulfilled (otherwise there appears a kinetic 
energy term). Two initial conditions cannot be simulta- 
neously fulfilled, if the equation contains time derivatives 
of the first-order only. Thus, there exists an initial stage 
of the dynamics of phase transitions (t < tmit)j which is 
not described by the standard Landau equation. 

For low velocities the heat transport is described by 
the heat conductivity equation cy = K AT, where k is 



2 



the heat conductivity and cy is the specific heat. Time 
scale of the temperature relaxation is tr — R 2 (tx)cy / k, 
where R(t) is the size of the germ. On the other hand, 
time scale of the density relaxation, following Eq. (Q]), is 
t p oc R (we show below that a germ of rather large size 
grows with constant velocity). Evolution of the germ 
is governed by the slowest mode. When sizes of germs 
begin to exceed the value i?f og , where i?f og is the size at 
which tx = t p , the growth is slown down. Thus number 
of germs with the size R ~ R[ og grows with time and 
there appears a metastable state called the fog. 

For further convenience we choose p v — p cr , T r = T cr 
and expand the Hclmholtz free energy in Sp and ST: 
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SF = F[p,T] - F[p T ,T t ]. Then a 
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less variables Sp = vip, 
we arrive at equation 

I = (2c/{\v 2 )Y'\ t = 2(dr, I + Q/(Xv 2 p r ), 
e = 2e/(Xv 3 ) = (8/3)\ST\^ [3 = cp 2 /[d Vl + G] 2 . 

Thus I oc \5T\- 1 ' 2 and t cx \ST\~ 1 . 

There exists an opinion, cf. Ref. [Tlj . that, if at some 
incident energy the trajectory passes in the vicinity of 
the critical point, the system may linger longer in this 
region due to divergence of susceptibilities that may re- 
flect on observables. Contrary, we argue that fluctua- 
tional effects in the vicinity of the critical point in heavy 
ion collisions can hardly be pronounced, since all relevant 
processes are proved to be frozen for ST — > 0, while the 
system passes this region during a finite time. 

To describe configurations of different symmetry we 
search two- phase solution of Eq. ^ in the form (J [2J , 

V> = Ttanh[£-£o(r)]+?/4, (4) 



£ = v£i + (2 + £3 for droplets/bubbles (d so i = 3), f = 
Vsi + £2 f° r r °ds (dsoi = 2) and £ = £1 = x/Z for kinks 
(dsoi = 1) in d = 3 space. For e > upper sign solution 
describes evolution of droplets (or rods and kinks of liquid 
phase) in a metastable super-cooled vapor medium. The 
lower sign solution circumscribes then bubbles (or kinks 
and rods of gas phase) in a stable liquid medium. 

The boundary layer has the length |£— £o( r )| ~ 1- Out- 
side this layer corrections to homogeneous solutions are 
exponentially small. Considering motion of the bound- 
ary for £o( T ) > 1 we may put £ ~ £0(1") m ©■ Then 
keeping only linear terms in e in Eq. ([3]) , we arrive at 
equation for £o( r ) : 
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Substituting Q in we obtain 

r rA , 27r^ 2 A 3 - d ^Xv i l d ^ 
SF[£ ] = 



r(4oi/2)r(i + (3 - 4oi)/2)p r 
x t?5o so 74oi + 2^- x /3 , 



(6) 



2A is the diameter, height of cylinder and the length of 
the squared plate for d so i = 3,2 and 1, respectively; T 
is the Euler T-function. The first term in ([6]) is the vol- 
ume term and the second one is the surface contribution, 
SF SUT {. At fixed volume in d = 3 space, the surface contri- 
bution for droplets/bubbles is smaller than for rods and 
slabs. Thereby if a germ prepared in a fluctuation is ini- 
tially nonspherical it acquires spherical form with passage 
of time. Surface term is SF snl . { = aS, S is the surface of 
the germ, a is the surface tension, and the gradient term 

S = h5F sur f. Thus we are 



in © is then 6F%£ 
able to find relations: a — ctq\ST\ 3 ^ 2 , a 2 = 32mn 2 r T cr c; 



31 p, 



I = 



There are two dimensionless parame- 
|): e and [3. The value ? distinguishes 



6T cr n cr |5T|!/2 

ters in Q and \ 
metastable and stable state minima in the free energy, 
[3 = (32T cr )^ 1 [di] T + Q]~ 2 cr 2 m controls dynamics. The 
larger viscosity and the smaller surface tension, the ef- 
fectively more viscous is the fluidity of germs. For [3 <C 1 
one deals with effectively viscous fluid and at f3 ^> 1, 
with perfect fluid. 

At hand of Eq. ([5]) consider analytically several typical 
solutions for germ evolution. Consider evolution of germs 
of stable phase in metastable matter. 

1) Short time evolution of a germ. For small r (initial 
stage) using Taylor expansion in r and assuming zero 
initial velocity, ^| T= o — 0, we obtain 

R(t) ~ Ro + (wt 2 /2) [1 - 2t/(3t [3)] 

valid for t « (^) V2 and * « t init = 2 Jk<±Ml oc 
Initial stage of the process proceeds with 



(dr ll +i I )\5T\ 

acceleration 



W = (d so l - 1)A« 2 (JZo - Rcr) /{RoRcr), 

which changes sign at the initial size Rq = R cr where 
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is the critical size. Germs with Rq < R cr shrink, 
while germs with Rq > R cr grow. For germs with 
\Rq — R cr \ *C R C r the size changes very slowly (w cx 
\5T\(Rq — R cr )/ R 2 r ). For undercritical germs of a small 
size, w oc — \5T\/Rq. Slabs of stable phase, being placed 
in a metastable medium, grow independently of what was 
their initial size. Note that the same value R cr follows 
from minimization of the free energy ©. 

2) Long time evolution of a large germ. For t ^> U n it, 
we may drop the term d 2 £,o/dr 2 in the l.h.s of Eq. ([5]). 
For R(t) ^> R cr , surface effects become unimportant and 
we arrive at the solution 



R(t) ~ R 



3\e\^/V2Xv i . 
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Germs grow with constant velocity. The time scale 
for the growth of the germ with size R ^S> R cr is 
t P = R/u &symp = (m/T CI .) 1 /2i?/((i8/?)i/2| ( 5 r |). Asymp- 
totic regime is reached at very large values of time, pro- 
vided the system is near the critical point. 

3) Long time evolution of a small germ. Describing 
germs of a small size (I <C R <C R cr , d so \ ^ 1) for t ^ 
^init, we can drop the term oc e in |5]). Then solution 
acquires the form 

R(t) ~ ^R 2 -2{d soX -l)tl 2 /t a . 

The time scale at which the initial germ of a small size 
dissolves is, t dis — 16n ^^( d ^+^) R o an j j s x phus, 

fluctuations of sufficiently small sizes are easily produced 
and dissolve rapidly. 

4) Fluctuations in spinodal region. Let the system 
be driven to a spinodal region where fluctuations of 
even infinitesimally small amplitudes and sizes may grow 
into a new phase. To demonstrate this we take the 
free energy SF to be close to its maximum (5F ~ 0). 
Then we linearize Eq. §3§ dropping ip 3 term. Setting 

-0 = — 4 + Re{7/>oe 7 ' , ' T+lfe ^} , ipo is an arbitrary but small 
real constant, we find two solutions, 



60MeV/fm 2 , Cmin = we evaluate /?gg> ~ 0.015g for 
(To — 40 MeV/fm 2 , that corresponds to the limit of effec- 



70(fc) = (-fc 2 ± ^ 4 + 8/3fc 2 - 4/3fc 4 )/(2/J). (7) 

Growing modes correspond to the choice of "+"-sign 
and k 2 < 2. The time scale at which an aerosol of 
germs develops is t acl — to/7^(fcm), k m corresponds to 
max{7^,(fc)}. For an effectively large viscosity (/3 <C 1) 
there are two solutions: the damped one, and the grow- 
ing one for k < \[2. The most rapidly growing mode is 
7V>(fcm) — 2, k m = 2/3 1 / 4 <C 1. The time scale charac- 
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4(di) +Co) 
9n cr T cr \8T\ ' 



tcrizing growth of this mode is 

The typical size of germs, R2 cr — Z / (2/? 1 / 4 ), increases 
with an increase of the viscosity. For k 2 > 2 both modes 
are damped. In the case of an effectively small viscos- 
ity (/? > 1) we get 7^(fc) ~ ±ki/2/l3y/l - k 2 /2, and 
7™ aa: (fcm = 1) = /3 -1 ^ 2 - The time scale characterizing 
growing modes, 4tr ~ h/li, = 2c 1/2 /(\v 2 ) oc ST' 1 , 
does not depend on the viscosity in this limit. The size 
scale of germs is R^ er ~ I. Modes with k 2 > 2 oscillate 
and do not grow into a stable phase. 

For the description of the hadron-QGP first-order 
phase transition we take values T cr ~ 162 MeV, 
n-cr/n^at — 1-3, as they follow from lattice calculations, 
see [13 ■ Parameters of the EoS are then as follows: 
8.76 • 10 2 (MeV • fm 3 ), b ~ 1.60 fm" 3 , A ~ 7.80 • 



1.56 • 10 4 g 2 |(5r|(MeV7fm b 

,3\ 



a ~ 

10- 5 q- 3 (fm7MeV ;! 

e ~ 2.02 • 10 2 (<5T) 2 (MeV/fm 3 ), where m is the effec- 
tive quark mass, q — (m/300MeV). Further we ob- 
tain l(T = 0) ~ 0.2 fm (radius of confinement) for 
cr ~ 40 MeV/fm 2 . If one used a ~ 100 MeV/fm 2 , 
one would estimate l(T — 0) ~ 0.5 fm. 

Next we use s ~ 7T 3 (T/T cr ) at T near T c p_cy ~ 
28T 3 (T/T cr ), as it follows from the lattice data [lj]. As- 
suming the minimal value of the viscosity 7y m i n = s /An ~ 



i 100 MeV/fni , 
0.2 < 1. Note 



tively very large viscosity. Even for Co 
to = 600 MeV we would get /?gg> - 
that following [l3[ the bulk viscosity diverges in the crit- 
ical point. If were so (/? — > 0), the quark- hadron system 
would behave as absolutely viscous fluid, like glass, in 



14j argues for a 



near critical region. Contrary, Ref. 
smooth behavior of the bulk viscosity. 

With (3 = 0.015 we further estimate t ~ 2| < 5T|- 1 fm, 
t p ~ 2.6i?g 1/2 |«T\ and i dis ~ 14gii (i? /fm). Typ- 
ical time for the formation of the aerosol is t^or — 
|(5T| _1 fm, and typical size of germs in aerosol is i?^ er ~ 
0.24|,5T|- 1 / 2 fm. Only i init ~ O.OS^TI- 1 fm proves 
to be small (excluding quite small 8T). Critical slow- 
ing down that limits growing of the a meson correlation 
length was discussed in [151 ]. 

For the thermal conductivity we use an estimation 
kqgp — aoij/m taking ao — 3 to recover the relation 
between values of k and rj for nuclear gas-liquid phase 
transition at low energies (l6| . The scale of the heat 
transport time is tx — 26q (i?/fm) 2 fm. Using that R cr ~ 
O.l^TI" 1 fm, we obtain R fog ~ Q.lq-^ 2 \8T\- X fm < 
R cr . The value i?f og proved to be very small (< 0.1 -f- 
1 fm). However, time scale tx is rather long. Therefore, 
the system most probably would have no time to fully de- 
velop a fog-like state in a hadron-quark phase transition 
in heavy ion collisions. 

For the system in the vicinity of the critical point all 
estimated time scales (except ij n it) are very large. If the 
system trajectory paths rather far from the critical point 
(T cr , Per), all time scales, except tx, become less than the 
typical life-time of the fireball (~ 10 fm at RHIC condi- 
tions). Reynolds numbers are < 1, being much smaller 
than the critical value (jj 1000). Thereby, turbulence 
regime is not reached. 

We solved numerically the general system of equations 
of nonideal hydrodynamics for d = 2. To illustrate the 
results we consider dynamics of overcritical and under- 
critical germs (disks) in infinite matter taking initial den- 
sity profile as p{x,y ; t = 0) = p out + (p in - Pout)&(Ra - 
r), r = y/x 2 + y 2 , pi n and p out are densities in stable 
and metastable homogeneous phases, respectively. 

In Fig. [1] we show the time evolution of a liquid disk 
(upper panel) and a gas disk (lower panel) for T/T cr — 
0.85. In the middle column we show dynamics of an ini- 
tially overcritical germ with R = 0.3L > R cr ~ 0.16L 
and in the right column, of undercritical germ Rq = 0.1L, 
L = 5 fm. The time snapshots are shown in Figure in 
units L. The configuration is computed for values of ki- 
netic parameters rj ~ 45 MeV/fm 2 and j3 ~ 0.2. We 
see that in case Rq > R cr (middle column) disks slowly 
grow with the time passage. For overcritical discs the ini- 
tially selected distribution acquires the tanh-like shape, 
see (UJ), only for t > (50 4- 100)L. Initial disks of a small 
size practically disappear for t ?S 10L = 50 fm. Due to 
the matter admission to the disk surface and the shape 
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FIG. 1: Isotherm for the pressure as function of the den- 
sity with initial and final states shown by dots (left column). 
Dash vertical line shows the Maxwell construction, MC. In 
the upper panel the initial state corresponds to the stable 
liquid phase disk in the metastable super-cooled gas, and in 
the lower panel, to a stable gas phase disk in the metastable 
super-heated liquid. Middle column shows time evolution of 
density profiles for the overcritical liquid disk (upper panel) 
and gas disk (lower panel). Numbers near curves (in L) are 
time snapshots. Right column, the same for initially under- 
critical liquid or gas disk. 

reconstruction, the density decreases in the liquid disk 
neighborhood below the value of the density in the ho- 
mogeneous metastable matter and it increases in the gas 
disk surrounding above the value of the density in the ho- 
mogeneous metastable matter (see the middle column). 

In Fig. [2] we demonstrate time evolution of the wave 
amplitudes, p(t) — p + Aof(t)sm(kr), for an undercritical 
value of the wave number k (left panel) and for an over- 
critical value (right panel) . In case of the overcritical k 



and effectively small viscosity (f3 = 20) we demonstrate 
change of the amplitude in the 3/2-periods of the oscilla- 
tion. Such a behavior fully agrees with that follows from 
our analytical treatment of the problem. 

Concluding, even in the spinodal region germs are 
growing slowly, if the system is somewhere in the vicinity 
of the critical point. Thus in heavy ion collisions the ex- 
panding fireball may linger in the QGP state, until T(t) 
decreases below the corresponding equilibrium value of 
the temperature of the phase transition. There exists 
a belief that strongly coupled QGP state, represents al- 
most perfect fluid [8[- We demonstrate the essential role 
of viscosity and surface tension in dynamics of first-order 
phase transitions, including the hadron-QGP one. Fluc- 
tuations in QGP (at a finite baryon density) grow and 
dissolve as if the fluid were very viscous. Variation of 
parameters in broad limits does not change conclusions. 




t/L t/L 



FIG. 2: Time evolution of wave amplitudes f(t) in aerosol 
for effectively small (/3 = 20, solid line) and large (/3 = 0.2, 
dash line) viscosity. Left panel: k — 21 /L (growing modes). 
Right: k = 81 /L (oscillation modes for large /3 and damping 
modes for small 0). Other parameters are the same as in Fig 
□ 
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